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Abstract. This paper studies algorithms similar to the Gaussian elimination algorithm in 
symplectic and split orthogonal groups. We discuss two applications of this algorithm in com¬ 
putational group theory. One computes the spinor norm and the other computes the double coset 
decomposition with respect to Siegel maximal parabolic subgroup. 


1. Introduction 


Gaussian elimination is a very old theme in Mathematies. It appeared in print as ehapter 
eight in a Chinese mathematieal text ealled, “The nine ehapters of the mathematieal art”. It is 
believed, a part of that book was written as early as 150 BCE. For a historieal perspeetive on 
Gaussian elimination, we refer to a niee work by Grear [j^. 

Due to many reasons 0. eomputational group theorists beeame interested in the construc¬ 
tive group recognition project. We will not go in any details of this projeet, but will refer an 
interested reader to the works of Leedham-Green and O’Brein [ |T3l and O’Brein [ [T^ Seetion 
9]. One ean also read a niee but slightly outdated review by Seress pO] , Matrix Groups] in this 
eontext. In dealing with eonstruetive group reeognition, one needs to solve the word problem 
in some generating set. As we know, in the speeial linear group SL( k), the word problem has 
an effieient solution in elementary transveetions - Gaussian elimination. 

In this paper, we work with Chevalley generators ||^ §11.3]. Chevalley generators for the 
speeial linear group are elementary transveetions. These Chevalley generators for other elas- 
sieal groups are known for a very long time. However, its use in row-eolumn operations in 
sympleetie and split orthogonal groups is new. We develop row-eolumn operations, very sim¬ 
ilar to the Gaussian elimination algorithm for speeial linear groups. We eall our algorithms 
Gaussian elimination in sympleetie and split orthogonal groups respeetively. Similar algorithm 


for twisted orthogonal groups and unitary groups p4| are being developed. 

The eurrent trend in computational group theory is to use standard generators. Using stan¬ 
dard generators, Brooksbank [|^, Costi Q solves the word problem in elassieal groups and 
Ambrose et. al. [[T| solves the membership problem in blaek-box groups. One advantage of 
using standard generators is that they are few in number. However, there is a disadvantages in 
working with them - they only work in finite fields. While working with Chevalley generators, 
our algorithms have two advantages: 

1. It works for arbitrary fields. 

2. It is mueh more effieient, as we demonstrate with an aetual implementation in Magma Q, 
see Figures 1 & 2. 
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From our algorithm, one can compute the spinor norm easily, see Section 5.2 Murray and 


Roney-Dougal [16| studied computing spinor norm earlier. Our algorithm can also be used 
to compute the double coset decomposition corresponding to the Siegel maximal parabolic 
subgroup, see Section 

Algorithms that we develop in this paper work only for a given bilinear form /3(see Equa¬ 
tions 2.2, 2.1 1 . Our algorithm work well on all characteristics for symplectic groups. However, 
for orthogonal groups our algorithms work only for odd characteristic. Henceforth, by suitable 
characteristics we mean all characteristics for the symplectic groups and zero or odd character¬ 
istic for orthogonal groups. It is known, in a field of suitable characteristics, all non-degenerate 
skew-symmetric bilinear forms are equivalent and all split (maximal Witt index) symmetric 
bilinear forms are equivalent. These equivalent bilinear forms are obtained by a base-change 
matrix and gives rise to conjugate linear groups. Though in our algorithm, we work with only 
one bilinear form /3, given by a fixed basis, with a suitable change of basis matrix our algorithm 
works on any equivalent bilinear forms. 

Another way to look at this paper, we have an algorithmic proof of this well-known theorem. 
For definitions of elementary matrices, one can look ahead to Section]^ 


Theorem A. Let k be afield of suitable characteristic. For d >A or I > 2 following holds: 

(£/) Every element of the split orthogonal group 0{d,k) can be written as a product of 
elementary matrices and a diagonal matrix. Furthermore, the diagonal matrix is of the 
form 

diag(l,..., 1, A, 1,..., 1, A G wheneverd = 2l 

diag(r?, 1 ,..., 1, A, 1,..., 1, A~^) XEk^andid = ±l whenever d = 21 + 1. 

Every element of the symplectic group Sp{2l,k) can be written as a product of elemen¬ 
tary matrices. 


This theorem has a surprising corollary. It follows that: 

Corollary B. In an orthogonal group, the image of X in k^ /k^^ is the spinor norm. 

This means that we have an efficient algorithm to compute the spinor norm. Since the com¬ 
mutator subgroup of the orthogonal group is the kernel of the spinor norm restricted to special 
orthogonal group, the above corollary is a membership test for the commutator subgroup in the 
orthogonal group. In other words, an element g in the special orthogonal group belongs to the 
commutator subgroup if and only if the A it produces in the Gaussian elimination algorithm is 
a square in the field, see Equation [2^ 

The bilinear form that we use and the generators that we define have its roots in the ab¬ 
stract root system of a semisimple Eie algebra and Chevalley groups defined by Chevalley and 
Steinberg [|6,21 1. However we assume no knowledge of Eie theory or Chevalley groups in this 
paper. 


2. Orthogonal and Symplectic Groups 

We begin with a brief introduction of orthogonal and symplectic groups. We follow Carter Q, 
Taylor [221 and Grove [ lOj in our introduction. In this section, we fix some notations which 
will be used throughout this paper. We denote the transpose of a matrix Xhy ^X. 

Eet y be a vector space of dimension d over a field k of suitable characteristic. Eet : V x 
V —)■ A be a bilinear form. By fixing a basis of V we can associate a matrix to /3. We shall abuse 
the notation slightly and denote the matrix of the bilinear form by /3 itself. Thus /3 (x, y) = 
^xfiy where x,y are column vectors. We will work with non-degenerate bilinear forms, which 
implies, det/3 7 ^ 0. A symmetric or skew-symmetric bilinear form j 8 satisfies /3 = ^/3 or /3 = 

2 












—^/3 respectively. In even order, for characteristic 2, the symmetric and skew-symmetric forms 
are the same. 


Definition 2.1 (Orthogonal Groups). A square matrix X of size d is called orthogonal if^X^X = 
/3 where j 8 is symmetric. The set of orthogonal matrices form the orthogonal group. 

In this paper, we deal with the split orthogonal group defined by one particular bilinear 
form defined over a field of zero or odd characteristics in Equations |2.2[ So any mention of 
orthogonal group means this one particular orthogonal group, unless stated otherwise. 

Definition 2.2 (Symplectic Group). A square matrix Xof size d is called symplectic if^XfiX = 
j8 where j8 is skew-symmetric. The set of symplectic matrices form the symplectic group. 


In this paper, we deal with the symplectic group defined by the bilinear form defined by 
Equation |2.1[ So any mention of symplectic group means this one particular symplectic group, 
unless stated otherwise. 

We write the dimension of V as J where d = 21 ox d = 21 and / > 1. In the case j 8 
is symmetric we define the corresponding quadratic form Q: V ^ khy Q{v) = jj5{v,v). Up 
to equivalence, there is an unique non-degenerate skew-symmetric bilinear form over a field 
k of suitable characteristics. Eurthermore a skew-symmetric bilinear form exists only in even 
dimension. We fix a basis of V as {ei,... so that the matrix j 8 is: 


( 2 . 1 ) 



The symplectic group with this /3 is denoted by Sp(2Z,k). 

Up to equivalence, there is a unique non-degenerate symmetric bilinear form of maximal 
Witt index over a field k of suitable characteristics. This is also called the split form. We 
fix a basis for odd dimension and {ei,.. . ... ,e_/} for even 

dimension so that the matrix /3 is: 


( 2 . 2 ) 


r 


P = { 


V 


/2 

0 

0 \ 



0 

0 


when d 

= 2 /+l 


II 

0 / 




0) 

1 

when d 

= 21. 


The orthogonal group corresponding to this form is a split orthogonal group. In this paper, 
we will simply call it the orthogonal group and this group will be denoted by 0{d^k). If we 
need to emphasize parity of the dimension, we will write 0(2/ -(- l,k) or 0(2/,k). We denote 
by kl{d,k) the commutator subgroup of the orthogonal group 0{d,k) which is equal to the 
commutator subgroup of SO{d,k). There is a well known exact sequence 


(2.3) 


0(J,k) 


SO(J,k) 



1 


m 

where 0 is the spinor norm. The spinor norm is defined as 0(g) = H Qi^i) where g = 

i=l 

Pvi '"Pvm is written as a product of reflections. Since the group SO(J,k) is of index 2 in 
0( k), we fix a generator for the quotient as wi = I — eij — 


3. Elementary Matrices and Elementary Operations 


In what follows, the scalar t varies over the field k and 1 <iJ< 1. Eurthermore, / > 2 which 
means d>A. We define teij as the matrix unit with t in the (/, j) position and zero everywhere 
else. We use eij to denote ^Cij. We often use the well known identity Cije^j = where 

dij is the Kronecker delta. 
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3.1. Elementary Matrices for 0(2/, k). We index rows by 1,2,— 1, —2,The 
elementary matrices are defined as follows: 


= 


for i ^ J, 

Xi,-j{t) = 


for i < J, 

= 


for i < J, 




and in matrix format 

0 ^ u r, , • / . 

Q r^-i j where R = l + tCij, i ^ j 
I R 


El: 


E2: ^ J where R is t{eij - cy/) for i < j 

E3: where R is t{eij - cy,) for i < j. 

The row and column operations: 

^R 0 
0 ^R-^ 

C D \0 ^R ^ 


ERl : 
ECl : 


A B 
C D 


( RA RB \ 

yR ^c t'r ^d) 

fAR B^R 
\CR D^'R ■ 


ERl: 


1 R 
0 7 


A B 
C D 


(a+rc b+rd\ 

\ C D ) 


ECl: 




(A AR + B\ 
l^C CR + DJ- 


ER3 : 



A B 
C D 


(A B \ 

\^7?A + C RB + DJ 


ECl: 



B 

D 


I 0 
R I 


/A + 57? B\ 
[c + DR DJ' 


3.2. Elementary Matrices for 0(2/ + l,k). We index rows by 0,1,...,/, — 1,..., —/. The 
elementary matrices are: 


Xi,j{t) = 


for / ^ 7 , 

Xi,-j{t) = 


for / < J, 

X-ijit) = 

1 

1 

+ 

for i < J, 

Xifi^t) = 



XQ,i{.t) = 

I + t{-le-ifi + eoj) - 



Written in matrix format, these four kind of elementary matrices are: 
/I 0 0 \ 

El: 10 7? 0 where 7? = 7 + tc,-^y; / 7 ^ j. 

\0 0 ^ 7 ? 7 
/I 0 0\ 

E2: 0 7 7? where 7? is t{ei j — ej i) ; i < j. 

\0 0 I 
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E3: 


E4a: 


E4b: 



0 where R is t{eij - ejj) ; i < j. 

V 

0 R \ 

I — ^7? j where R = tei 

0 / J 

R 0\ 

I 0 where R = tei 
-^RR l) 


Here ei is the row veetor with 1 at place and zero elsewhere. 


3.3. Elementary operations for 0(2/ + l,k). Eet g = 


a X 
E A 
F C 


Y 

B 

D 


bea (2/ + 1) x (2/+1) 


matrix where A, 5, C,D are I x I matrices. The matrices X = (Xi,X 2 ,... ,X/), 7 = (71,72, • • • j^/), 
E = ^(£'i,£' 2 , • • ■ ,E/) and F = ^(Fi,F 2 ,... ,F/). Eet a G k. Eet us note the effect of multiplica¬ 
tion by elementary matrices from above. 


fl 0 0 \ fa X Y\ fa X Y \ 

ERl: lo R 0 E A 5 = RA E5 

\0 0 7 \F c dJ yR ^F ^R ^r ^dJ 


ECl : 



/I 0 
0 R 
\0 0 



fa XR Y^R y 
E AR B^R ^ 
\F CR 7 


fl 0 0\ fa X Y\ fa X Y \ 

ERl : 0/E EA5 = E+RE A + RC B + RD 

\0 0 ij \F C DJ \ E C D ) 


fa X Y\ fl 0 0\ 
ECl: \E a 5 0 1 5 

\F C dJ \0 0 ij 


fa X XR + Y\ 
E A AR + B 
\F C cr + dJ 


fl 0 0\ fa X Y\ 
ER3 : 0 I 0 E A 5 

\0 R IJ \F C d) 

fa X Y\ fl 0 0\ 
EC3: E A 5 0 I 0 

\E C dJ \0 R l) 


fa X 7 \ 

E A 5 

\EE-fE EA + C RB + dJ 

fa X + YR 7\ 

E A+BR B . 

\F c + dr d) 


Eor E4 we only write the equations that we need later. 

• Eet the matrix g has C = diag( 

ERA : [(/ + teo_i)g]o,i = Xi A- tdi 


ECA : [g{I + teo,-/ - 2te/,o - _,)] = Fi- 2tdi. 

• Eet the matrix g has A = diag(,..., J/). 

EE4 : [(/ -h teo^i - Ite^i^o - t^e-ij)g]o^i = Xi + tdi 

ECA : [g{I + teo-i - 2te/,o - t^ei~i)]ifi = Ei - ltdi. 
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3.4. Elementary Matrices for Sp(2/,fc). 


Xij{t) = 


for i ^ 7, 

= 

I + t{ei-j + ej-i) 

for i < 7, 

X-i,j{t) = 

k + t{e-ij + e^j^i) 

for i < 7, 

= 

I + tei-i 



13 - 



There are three kinds of elementary matriees. 

El: wherei?=/ + tc,j;iV7. 

E2: where R is either t{eij + ej^i); i < j or tei^i. 

E3: f ^ ^ j where R is either t{eij + cy;); i < j or tei^i. 


3.5. Elementary Operations for Sp(2/,k). Eet ^ be a 2/ x 21 matrix written in 

block form of size / x Z. Then the row and column operations are as follows: 


ER\ : 




A B 
C D 


RA 

^R~^C 


RB \ 
'^R ^D) 


ECl : 


A B 
C D 



fAR B'^R-^\ 

d'^r-^) ■ 


ER2 : 


1 R\ A B 


0 1 \C D 


(A+RC B + RD\ 
\ C D ) 


ECl: 




fA AR + B\ 
[c CR + Dj- 


ER3 : 



0\ [A B 
I VC D 


A B \ 

RA + C RB + DJ 


EC3 : 



B 

D 


I 0 
R I 


fA + BR B\ 
[c + DR DJ- 


Remark: We would require row-interchange of Z* row with —Z* row in our algorithms. In 
the case of symplectic and odd-order orthogonal groups this row-interchange is a product of 
elementary matrices. In the case of even-order orthogonal groups, we need to add to generators 
a row-interchange matrix w/. Eor more see Eemma[43j 


4. Gaussian Elimination in orthogonal and symplectic groups 

Recall the field k is of suitable characteristic. Cohen, Murray and Taylor [[7| proposed a 
generalized row-column operations, using a representation of Chevalley groups. The key idea 
there was to bring down an element to a maximal parabolic subgroup and repeat the process 
inductively. The emphasis there was to represent generators as symbols so that it takes less 
memory to store. Here we use the natural matrix representation of these groups. 
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4.1. Gaussian Elimination for Sp(2Z, A:) and 0(2Z,A:). The algorithm is as follows. 


Step 1: Input: A matrix g = 


C 


Output: The matrix gi = 


D 

Ai 

Cl 


whieh belongs to Sp(2Z,fc) or 0(2Z,fc). 
is one of the following kind: 


Di 


a: The matrix A i is a diagonal matrix diag(l,..., 1, A) with A 7 ^ 0 andCi is 


to m and Ci is of the form 


where Cn is an m x m symmetrie matrix 


Cii C 12 
C 21 C 22 

where Cn is symmetrie when g is in Sp(2Z,A) and skew-symmetrie when g is in 
0(2Z, A) and is of size Z — 1. Furthermore, C 12 = A^C 2 i when g is in Sp(2Z,A) and 
C 12 = —A^C 2 i, C 22 = 0 when g is in 0(2Z,A). 

The matrix Ai is a diagonal matrix diag( 1,..., 1,0,..., 0) with number of Is equal 

Cii 0 
C 21 C 22 

when g is in Sp(2Z, A) and skew-symmetrie when g is in 0(2Z, A). 

Justification: Observe the effeet of ERl and ECl on the bloek A. This amounts to 
Gaussian elimination on a Z x Z matrix. Thus we ean reduee A to a diagonal matrix and 
Corollary |4.2| makes sure that C has required form 

Ai Bi' 

Cl Di 

M B 2 

0 ^A 2 

Justification: Observe the effeet of ER3. It ehanges Ci to RAi -l-Ci. Using Eemma [43] 
we ean make the matrix Ci the zero matrix in the first ease and Cn the zero matrix in 
the seeond ease. Eurthermore, in the seeond ease, we use Eemma [A 6 ] to interehange the 
rows so that we get the zero matrix in th e plaee of Cn If required use ERl and ECl to 

4.4 ensures thatDi beeomes ^AA^ 


Step 2: Input: matrix gi = 
Output: matrix g 2 = 


; A 2 is a diagonal matrix diag( 1 ,..., 1 , A). 


make Ai a diagonal matrix. Eemma 
^A2 B 2 
0 ^A- 


Step 3: Input: matrix g 2 = 


Output: Matrix g^ = 


A2 

0 


0 


; A 2 is a diagonal matrix diag( 1 ,..., 1 , A). 
); A 2 is diagonal matrix diag( 1 ,..., 1 , A). 


Justification: Using Corollary 4.3 we see that the matrix B 2 has eertain form. We ean 
use ER2 to make the matrix B 2 a zero matrix beeause of Eemma |4.5[ 

The algorithm terminates here for 0(2Z, A). However for Sp(2Z, A) there is one more step. 

Step 4: Input: matrix g^ = diag(l,...,1,A,A^^). 

Output: Identity matrix 

Justification: This can be written as a product of elementary matrices by the first part 
of Eemma 14.71 

4.2. Gaussian Elimination for 0(2Z -I- 1, A). An overview of the algorithm is as follows: 

/a X y' 

Step 1: Input: matrix g= \e A B which belongs to 0(2Z-|-1, A); 

\F C 


of one of the following kind: 



a 

Output: matrix gi = | Ei 

a: The matrix A 1 is a diagonal matrix diag(1,..., 1, A) with A 7 ^ 0. 
b: The matrix Ai is a diagonal matrix diag( 1,..., 1,0,..., 0) with number of Is equal 
to m and m < Z. 
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Justification: Using ERl and ECl we do the elassieal Gaussian elimination on a / x Z 
matrix A. 

/a Xi Yi 

Step 2: Input: matrix gi = Ei Ai Bi 

\Fi Cl Di 
(a2 X 2 Y 2 

Output: matrix g 2 = E 2 ^2 B 2 | of one of the following kind: 

\F2 C 2 E>2^ 

a: The matrix A 2 is diag(l, 1,..., l,-^) with A 7 ^ 0, X 2 = 0 = E 2 and C 2 is of the form 

where Cn is skew-symmetrie of size / — 1. 

b: The matrix A 2 is diag( 1,..., 1,0,..., 0) with number of Is equal to m; X 2 and E 2 

Cn 0 


have first m entries 0, and C2 is of the form 


C21 C22 


where Cn is an m x m 


skew-symmetrie. 

Justification: Once we have Ai in diagonal form we use ER 4 and EC 4 to change Xi 
and El in the required form. Then Eemma | 4 . 8 1 makes sure that Ci has required form. 
[0C2 X2 Y2 

Step 3 : Input: matrix g2= E2 A2 B2 

\F2 C 2 D2j 

Output: 


a: matrix g3 = 


b: matrix g3 


a3 0 73 

0 A3 E3 
^Ft, 0 

«3 X3 73' 
E3 A3 E3 
E3 C3 D3 


where A3 is diag( 1A). 


where A3 is diag( 1,..., 1,0,..., 0) with number of Is 


equal to m; ^3 and E3 have first m entries 0 , and C3 is of the form 


with A4 diagonal matrix diag( 1,..., 1, A). 


0 0 
^C 2 l € 22 ^ 

Justification: Observe the effect of ER3 and the Eemma |43| ensures the required form. 
/a3 X 3 73 ' 

Step 4: Input: g 3 = E 3 A 3 E 3 

\E 3 C 3 

/t} 0 0 

Output: g 4 = j 0 A 4 B 4 
\0 0 A 4 I 

Justification: In the first case, Eemma |4. 10| ensures the result. In the second case we 
interchange Z* with — Z* for m + 1 < i < 1. This will make C 3 = 0. Then if needed we 
use ERl and ECl on A 3 to make it a diagonal. Eemma [439] ensures that A 3 has full rank. 
Eurthermore, we can use ER4 and EC4 to make X 3 = 0 and E 3 = 0. Eemma 4.10 gives 
the required form. 

/i} 0 0 \ 

Step 5: Input: g 4 = j 0 A 4 E 4 j withA 4 = diag(l,..., 1, A). 

VO 0 A^V 

Output: g 5 = diag(r^ 1... 1,A,A^y. 

Justification: Eemma |4T^ ensures that B 4 is of certain kind. Use ER2 to make B 4 = 0. 
Eemma [43] ensures that the first diagonal entry is ■&. 
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Proof of the Theorem A. {s^) Let g E 0{d,k). Using the above algorithm we can reduce g 

to a diagonal matrix of the required form. 

If g G Sp(2/,A:) using above algorithm and Lemma 4.7 we can reduce g to the identity. 
Thus g is a product of elementary matrices. 


This algorithm can be used to compute the determinant of an orthogonal group. In the even- 
order case, refer to Algorithm 4.1. After Step 1, either Ai has full rank, or the rank of Ai is 
m. In the first case the determinant of g is 1 and in the second case it is (—1)^^'". It is fairly 
straightforward to see this. Observe that all elementary matrix other than the row-interchange 
matrix has determinant one and the number of row-interchange that we do is / — m, where m is 
the rank of Ai. 

In the odd-order case, as we have observed that the row-interchange matrix is a product of 
elementary matrices and all elementary matrices have determinant one. It then follows clearly 
that the determinant of g is r?. 


4.3. Time-complexity of the above algorithm. We establish that the worst case time-complexity 
of the above algorithm is 0 (/^). 

In Step I, we make A a diagonal matrix by row-column operations. That has complexity 

0(P). 

In Step 2, Cl -|-7?Ai is multiplying two rows by a field element and two additions. In the 
worst case, it has to be done 0(1) times and done 0(1^) many times. So the complexity 
isO(/3). 

Step 3 is similar to Step 2 above and has complexity 0(/^). 

Step 4 has only a few steps that is independent of 1. 

Then clearly, the time-complexity of our algorithm is 0(/^). 

4.4. Lemmas used in the justification of the Gaussian elimination. To justify the steps of 
Gaussian algorithm we need several lemmas. Some of these might be well known to experts 
but we include them here for the convenience of the reader. 


Lemma 4.1. Let Y = diag(l,..., 1, A,..., A) be of size I with number of U equal to m < 1. 
Let X be a matrix of size 21 such that YX is symmetric (skew-symmetric) then X is of the form 

/All -^^21 I Xii is symmetric (skew symmetric) and A12 = A^2i (2(^X2 = — 

V^ 2 i A22 J 

Furthermore, ifXf^O thenX 22 is symmetric (skew-symmetric). 

Proof. We observe that the matrix YX = ^ • The condition that YX is symmetric 

implies An (and A 22 if A 7 ^ 0) is symmetric and A 12 = A^A 2 i. • 

A B 


Corollary 4.2. Let d] either in Sp{2l,k) or 0(2/, k). 


(1) If A is a diagonal matrix diag(l,..., 1,0,... ,0) with number of Is equal to m(< /) then 

'Cii 0 

fd2l C22^ 

skew-symmetric if g is orthogonal. 


the matrix C is 


where Cn is an mx m symmetric if g is symplectic and is 


(2) If A is a diagonal matrix diag(l, 1,..., 1, A) then the matrix C is 


Cn A'C2i 


where 


fd2i C22 

Cii is an (/ — 1 ) x (/ — 1 ) symmetric if g is symplectic and Cn is skew-symmetric with 
C 22 = 0 ifg is orthogonal. 
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Proof. We use the condition that g satisfies = /3 and get AC is symmetric (using A = ^A as 


A is diagonal) when g is symplectic and skew-symmetric when g is orthogonal. The Lemma 4.1 
gives the required form for C. • 


Corollary 4.3. Let g = 


where A = diag(l ,... ,l,X)be an element of either Sp{2fk) 


or 0(21,k) then the matrix B is of the form ( 


where B\ \ is a symmetric matrix 


^A B 
0 A-\ 

B 22 

of size I —I if g is symplectic and is skew-symmetric along with B 22 = 0 ifg is orthogonal. 

Proof. Yet again, we use the condition that g satisfies = /3 and A = ^ to get A^^B is 
symmetric if g is symplectic and is skew-symmetric if g is orthogonal. Then Lemma [4~T] gives 
the required form for 5. • 


Lemma 4.4. Let g = EGL(2fk). Then, 

(1) the element g belongs to Sp(2l, k) if and only ifD = ^A^^ and BA = A^B and 

(2) the element g belongs to 0(21, k) if and only ifD = ^A^^ and BA — —A^B. 

Proof. This follows by simple computation using = /3. • 


Lemma 4.5. Let Y = diag(l, 1,..., 1, A) be of size I where A 7 ^ 0 and X = (xij) be a matrix 
such that YX is symmetric (skew-symmetric). Then X = (i?i -f + ■ • 0^ where each Rm is of 
the form t(eij -1- ejf for some i < j or of the form tetjfor some i (in the case of skew-symmetric 
each Rffi is of the form t(eij — ej^i)for some i <j)- 


Proof. Since YX is symmetric, the matrix Y is of the form 


^11 

^21 


^12 

xii 


where Yii is symmetric 


and Y 21 is a row of size / — 1 and Y 12 = A ^ 21 - Clearly any such matrix is sum of the matrices 
of the form RmY■ A similar calculation proves the result for the skew-symmetric case. • 


We need certain Weyl group elements which can be used for interchanging rows. We use a 
formula Wr = x^(l)x_r( —l)xr(l) from the theory of Chevalley groups 0 Lemma 6.4.4] and 
construct elements Wij and Wi-j. In our algorithm, we need elements which interchanges i* 
row with —i* row for any i. The element wt-t that we create when multiplied to a matrix g 
interchanges its rows while simultaneously multiplying some rows by — 1. However that does 
no harm to our algorithm. 


Lemma 4.6. For I <i <1, 

(1) the element wi-i = I-\- — ei^i — e^i-i G Sp(2l, k) is a product of elementary 

matrices. 

(2) the element wi-i = I — 2eo,o ~ ~ ~ — ^-i,i ^ 0(21 -|- 1, A) is a product of 

elementary matrices. 

(3) The element Wi-i = I — et-i — — e^i-i G 0(21, k) is a product of elementary 

matrices. 


Proof. For the symplectic group Sp(2Z,A) we have w, -,- = Xi-i(l)x-ij(—l)xi-i(l). For the 
orthogonal group 0(2/ -b 1, A) we have w/= xo,,(-l)^i,o(l)j‘^o,i(-l)- 

For the orthogonal group 0(21,k) we inductively produce these elements. First we get 


W; 




j — (1 + ei-j -ej-i)(l + e^ 


ij 


:j)(I-\-ei-j ^i,i 


' ^jJ 


■ep-i + e^ 




■e^j,iandwij = I- 




-e^i^^i-e- 


-j-j + 


Now we set w/^_/ = w;. Then compute 

^-(Z-l),-(/-l) - ^(Z-l),-(/-l) - e-(Z-l),(Z-l)- 
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Lemma 4.7. (1) In the case of Sp{2l,k), the element diag(l,..., lA, 1,..., 1, is a 

product of elementary matrices. 

(2) In the case ofO{2l + 1, A:) the diagonal element diag( —1,1, A^^) is a 
product of elementary matrices. 

Proof. In the case of Sp(2/,fc), we compute = (/ + te;__/)(/ —+ = 

^ — ^i.i — ^-i,-i + tei-i—t~^e^ij and then compute /z/(A) = w/^_/(A)w/__/(—1) which is the 
required element. 

In the case of 0(2/ + 1,A:) we compute wip{t) = xoj{—t)xip{t~^)xoj{—t) =7 — — 

t^e-ij — ei^i — 2eo,o — and multiply it with w/^o(l) to got the required matrix. • 


Remark : In the case of 0(21,k), the element diag(l,..., lA^, 1,..., 1, A and in the case of 

0 ( 21 +I,k) the element diag( 1 , 1 ,..., 1 A^,A is a product of elementary matrices. 


/a X 

Lemma 4.8. Let g = j * A * be in 0{2l-\-l,k). 
\ * C * I 


(1) If A = diag(l,..., Ij X = 0 then C is of the form 
skew-symmetric. 

(2) If A = diag( 1,..., 1,0,..., 0) with number of li equal m < I and X has first m entries 0 

then C is of the form J with Cn skew-symmetric. 


iJCA + ^AC). In the first case X = 0, 
4.2| to get the required form for C. In the second case, we note that ^XX has top-left 


Proof. We use the equation ^gfig = fi and get 2^XX = 
so use 


block 0 and get the required form. 


/a X Y\ 

Lemma 4.9. Let g = j * A *1 be in 0{2l-\-l,k) then X = 0 and D = ^A^^. 

0 dJ 


Proof. We compute ^gfig = /3 and get 2^XX = 0 and 2^XY -|- ^AD = 7. This gives the required 
result. • 

fa 0 Y\ 

Lemma 4.10. Let g = I 0 A B \, with A an invertible diagonal matrix. Then, g G 0(21 -|- 

\F 0 dJ 

\,k) if and only if Ot} = \,F = 0 = Y, D =A~^ and ^DB + ^BD = 0. 


Proof. 




fa 0 ^F\ f2 0 0\ 

0 ^A 0 0 0 7 

yY ^B ^dJ \0 I oj 



2 a^ 

^AF 

af'Y + ^BF 


t'FA 

0 

^DA 


2aY + ^FB '' 
^AD 

pryy I I Tdi^ 


Equating this with /3 gives us the required result. 
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5. Computing Spinor Norm for orthogonal groups 


In this section, we show how we can use Gaussian elimination to compute spinor norm 
for orthogonal groups. The classical way to define spinor norm is via Clifford algebras pO} 
Chapters 8 & 9]. Spinor norm is a group homomorphism 0: 0{d,k) —)■ /k^^, restriction 
of which to SO{d,k) gives Q.{d,k) as kernel. However, in practice, it is difficult to use that 
definition to compute the spinor norm. Wall [231, Zassenhaus and Hahn [ 111 developed a 
theory to compute the spinor norm. For our exposition, we follow [ [22l Chapter 11]. 

Let g be an element of the orthogonal group. Consider g as a linear transformation. Futher- 
more, denote g = I — g and Vg= g{V) and = ker{g). Using /3 we define Wall’s bilinear form 
[, ]g on Vg as follows: 


[u,v]g = P{u,y), where, v = g{y). 
This bilinear form satisfies following properties: 

(1) [M,v]g + [v,M]g = /3 (m,v) and [u,u\g = Q{u) for all m,v G Vg. 

(2) g is an isometry on Vg with respect to [, ]g. 

(3) [v, u\g = — [m, gv] for all u,v &Vg. 

(4) [, ]g is non-degenerate. 


Then the spinor norm is 

0 (g) = disc(yg, [, ]g)ifg^/ 

extended to I by defining 0(7) = 1. An element g is called regular if Vg is non-degenerate 
subspace of V with respect to the form /3. Hahn [11, Proposition 2.1] proved that for a regular 
element g the spinor norm is 0(g) = det(g|u )disc(yg). This gives. 


Proposition 5.1. (1) Fora reflection Pv, 0(Pv) = Q{v). 

(2) 0(-l) = JAc(y,/3). 

(3) For a unipotent element g the spinor norm is trivial, i.e., 0(g) = 1. 


Murray and Roney-Dougal [ 16 j used the formula of Hahn to compute spinor norm. However, 
we show (Corollary B) that the Gaussian elimination developed in Section [^ouptuts the spinor 
norm quickly. First we observe the following: 


Lemma 5.2. For the group 0{d,k), d >A, 

(1) &{xij{t)) = ®{x-ij{t)) = &{xi-j{t)) = 1. Furthermore, in odd case we also have 

®{xi,o{t))_= T = 0(ji;o,,(O)- 

(2) 0(w/) = 1. 

(3) 0 (diag(l,..., 1 ,A,A^^)) = A. 

Proof. We use Proposition [5.1| The first claim follows from the fact that all elementary matrices 
are unipotent. The element w/ = is a reflection thus 0(w/) = Q{ei -|- e_/) = 1. 

For the third part we note that diag(l,..., 1, A, 1,..., 1, A^^) = p(e,+e_,)p(e,+Ae_,) and hence 

thespinornorm0(diag(l,...,l,A,l,...,l,A^i)) = 0(p(g,+;Le_,)) = + =^- • 

Proof of Corollary B. Let g G 0{d, k) . From Theorem A, we write g as a product of elementary 
matrices and a diagonal matrix diag( 1,...,1,A,A^^). Furthermore, the spinor norm of 
an elementary matrix is 1 . Thus 0(g) = 0 (diag(l,..., 1, A, 1 ,..., 1 , A^^)) = A. • 


6 . Double coset decomposition eor Siegel maximal parabolic 

In this section, we compute the double coset decomposition with respect to Siegel maximal 
parabolic subgroup using our algorithm. Let P be the Siegel maximal parabolic of G where 
G is either 0{d,k) or Sp(2Z,A). In Lie theory, a parabolic is obtained by fixing a subset of 
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simple roots 0 Section 8.3]. Siegel maximal parabolic corresponds to the subset consisting of 
all but the last simple root. Geometrically, a parabolic subgroup is obtained as fixed subgroup 
of a totally isotropic flag p3| Proposition 12.13]. The Siegel maximal parabolic is the fixed 
subgroup of following isotropic flag (with the basis in Section |^: 

{0}c{ei,...,ez}cy. 

a 0 Y" 

E A B 

/ 0 D, 

The problem is to get the double coset decomposition P\G/P. That is, we want to write 
G = IJ P(0P as disjoint union where IT is a finite subset of G. Equivalently, given g G G v^e 
coeW 

need an algorithm to determine the unique CO GW such that g G PcoP. If G is connected with 
Weyl group IT and suppose Wp is the Weyl group corresponding to P then [|^ Proposition 2.8.1] 

P\G/P <—^ Wp\W/Wp. 


Thus P is of the form 


in 0 ( 2 / + l,fc) and 


B 

D 


in Sp(2/,/:) and 0(2/,/:). 


We need a slight variation of this as the orthogonal group is not connected. 

In the case of Sp(2/,A:), the Weyl group IT =N{T)/T where T is a diagonal maximal torus 
and T = {diag(Ai,..., A/, Af \. • •; ^The group W is isomorphic to a subgroup 

of S 21 , the symmetric group on 2 / symbols {1 ,..., /, — 1,..., —/} and is generated by elements 
and wz,-,-which map to permutations (/,/+1)(—/,—(/+1)) and (/,—/) respectively. Thus 
W is isomorphic to Si xi (Z/2Z)^ and the subgroup Wp is generated by wz,z+i which proves that 
the subgroup {(/,/+ !)(—/,—(/+ 1)) | 1 < / < /} is isomorphic to 5/. 

For Sp(2/,A), we set IT = {(Oq = /,tOz = wi _i ■■■Wi-i | 1 < / < /} and note that IT = 

U Wp(oW/p. 

i =0 

In the case of 0(J, A), we set IT = {(Oo = £0/ = wi-i ■ ■■Wi-j | 1 < / < /} where w/is 

inductively produced (see Lemma [43] ). 

Theorem 6.1. Let P be the Siegel maximal parabolic subgroup in G, where G is either 0{d, k) 
or Sp{d, k). Let g G G. Then there is an efficient algorithm to determine CO such that g G PooP. 
Furthermore, IT the set of all 00 s is a finite set of 1+I elements where d = 21 or 21+1. 


Proof In this proof we proceed with a similar but slightly different Gaussian elimination algo- 

B\ 


rithm. Recall that g = 


A 

C 


DJ 


whenever g belongs to Sp (2/, A) or O (2/, A) or g = | E A 

.F C 


there are elements p\,P 2 ^P such that pigp 2 = 


of A\ zero, i.e., there is a p 3 G P such that P 3 Pigp 2 = 


whenever g belongs to 0(2/ 1, A). In our algorithm, we made A into a diagonal matrix. In¬ 

stead of that, we can use elementary matrices ERl and ECl to make C into a diagonal matrix 
and then do the row interchange to make A into a diagonal matrix and C a zero matrix. If we 
do that, we note that elementary matrices El and E2 are in P. The proof is just keeping track 
of elements of P in this Gaussian elimination algorithm. The step 1 in the algorithm says that 

Bi 

,Ci 

m non-zero entries. Clearly m = 0 if and only if g G P. In that case g is in the double coset 
PcOqP = P. Now suppose m > 1. Then in Step 2 we multiply by E2 to make the first m rows 

Pi' 

,Ci Di^ 

are zero. After this we interchange rows i with —/ for 1 < / < m which makes Ci zero, i.e., 

''M P2 
0 P> 2 , 
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where Ci is a diagonal matrix with 


where first m rows of Ai 


multiplying by we get cOmPiPigpi = 


G P. Thus g G PoOmP. 





For 0(2/ + 1, A:) we note that the elementary matriees El, E2 and E4a are in P. Rest of the 
proof is similar to the earlier ease and follows by earefully keeping traek of elementary matrices 
used in our algorithm in Section [43} • 


7. Conclusions 


We conclude this paper with some implementation results of our algorithm and some com¬ 
parisons with an existing algorithm. Before we state those, we wander a little towards a di¬ 
rection for further research. Chevalley generators are known for a long time. However, its use 
in row-column operations is new. Earlier, in computational group theory, while working with 
quasisimple groups, the generators of choice were always the standard generators. This use 
of Chevalley generaotrs can bring in a paradigm shift with algorithms in matrix groups. It is 


now a very interesting project to redo the constructive group recognition project [ 131 with our 
algorithms and Chevalley generators. 

Now some implementation results, we implemented our algorithm in magma Q. We found 
our implementation to be fast and stable. In magma, Costi and C. Schneider installed a function 
ClassicalRewriteNatural. It is the row-column operation developed by Costi @ in natural 
representation. We tested the time taken by our algorithm and the one taken by the Magma 
function. To do this test, we followed Costi ||^ Table 6.1] as closely as possible. Two kind of 
simulations were done. In one case, we fixed the size of the field at 7^® and varied the size of 
the matrix from 20 to 60. To time both these algorithms for any particular input, we took one 
thousand random samples from the group and run the algorithm for each one of them. Then the 
final time was the average of this one thousand random repetitions. The times were tabulated 
and presented below. 

In the other case, we kept the size of the matrix fixed at 20 and we varied the size of the field, 
keeping the characteristic fixed at 7. In many cases the magma computation for the function 
ClassicalRewriteNatural will not stop in a reasonable amount of time or will give an error and 
not finish computing. In those cases, though our algorithm worked perfectly, we were unable 
to get adequate data to plot and are represented by gaps in the graph drawn. Here also the times 
are the average of one thousand random repetitions. It seems that our algorithms perform 


Size of matrix 20 Size of the field 7 




Eigure 1. Some simulations comparing our algorithm with the one inbuilt in 
Magma for even-order orthogonal groups 


better than that of Costi’s on all fronts. 
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Size of the matrix 20 


Size of the field 7 
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Figure 2. Some simulations comparing our algorithm with the one inbuilt in 
Magma for symplectic groups 
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